Threedimensional Local Porosity Analysis of Porous Media 



B. Biswal 1 ' 2 , C. Manwart 1 and R. Hilfer 1 ' 3 

1 ICA-1, Universitat Stuttgart, Pfaffenwaldring 27, 70569 Stuttgart 

2 Department of Physics & Electronics, Sri Venkateswara College, 

University of Delhi, New Delhi - 110 021, India 

3 Institut fur Physik, Universitat Mainz, 55099 Mainz, Germany 

A quantitative comparison of the pore space geometry for three natural sandstones is presented. 
The comparison is based on local porosity theory which provides a geometric characterization of 
stochastic microstructures. The characterization focusses on porosity and connectivity fluctuations. 
Porosity fluctuations are measured using local porosity distributions while connectivity fluctuations 
are measured using local percolation probabilities. We report the first measurement of local percola- 
tion probability functions for experimentally obtained threedimensional pore space reconstructions. 
Our results suggest to use local porosity distributions and percolation probabilities as a quantitative 
method to compare microstructures between models and experiment. 
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I. INTRODUCTION 



A large number of microscopic models have been proposed to represent the microstructure of porous media [jl] 12| . 
Representative microscopic models are a prerequisite for studying transport properties such as fluid flow or sound 
propagation in oil reservoirs, aquifers or other random media. 

Microscopic models are not unique, and hence it is necessary to have criteria for comparing them among each other 
and to the experimental porous microstructure jl3] [l?]] . This is particularly important for attempts to generate porous 
microstructures in an automatic computerized process [^8|j7|J^,^9| , or to decide quantitatively whether the connectivity 
is percolationlike as is often assumed in models |2fj| , pT| , p| . 

Detailed microscopic models contain so many geometrical features that the rough comparison based only on porosity 
and specific surface is insufficient. The problem is to find general geometric characterization methods to test how well 
a model represents the microstructure found in reality. Given such tools they can then be used to constrain the input 
parameters of the models. 

General geometric characterization methods traditionally employ only porosities, specific surface areas, and some- 
times correlation functions ||,|22|,^3],f7 24 1 . Recently novel tools based on local porosity theory became available for 



the comparison of stochastic microstructures p5|- f^fi"3|]2"8| |3(],|l5| . Local porosity theory is currently the most general 
geometric characterization method because it contains as a special case also the characterization through correlation 
functions (see [[l5) for details). 

Local porosity theory contains two geometric characteristics. The first are local porosity distributions, the second 
are local percolation probabilities |l5f| . While local porosity distributions have been measured previously on artificial 
and real samples p^j29f -|3l|] no reliable measurement has been made up to now for the local percolation probabilities. 
The main impediment has been the absence of accurate threedimensional pore space representations for real rocks. 

The objective of the work reported here has been to measure the local percolation probabilities of natural sandstones 
thereby providing a new geometric characteristic against which microscopic pore space models can be compared. 



II. MEASURED QUANTITIES 
A. Local Porosity Distributions 



Local porosity distributions were originally introduced as a quantitative substitute for pore size distributions [|25| . 
The idea is to measure porosity or other well defined geometric observables within a bounded (compact) subset of the 
porous medium and to collect these measurements into various histograms (empirical probability densities). 

Imagine a porous medium occupying a subset § C M. d of the physical space (d — 3 in the following). For the 
data analysed here the set § is a rectangular parallelepiped whose sidelengths are Mi , Mi and M3 in units of the 
lattice constant a (resolution) of a simple cubic lattice. The sample S contains two disjoint subsets § = P U M with 
P n M = where P is the pore space and M is the rock or mineral matrix and is the empty set. In practice the 
sample is discrctized, and the configuration of the two sets P and M is given as an Mi x M2 x M3-array of two numbers 
representing P and M respectively. Let K(x, L) denote a cube of sidelength L centered at the lattice vector x. The 
set K(x, L) defines a measurement cell inside of which local geometric properties such as porosity or specific internal 
surface are measured. The local porosity in this measurement cell K(x, L) is defined as 

VV ' ' V(K{x,L)) V ; 
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where V(<G) is the volume of the set G C M. d . The local porosity distribution fi((f>, L) is defined as 

^,£) = iV^-^I)) (2.2) 

TO * — ' 
x 

where m is the number of placements of the measurement cell K(x, L) . The results presented below are obtained 
by placing K(x, L) on all lattice sites x which are at least a distance L/2 from the boundary of S, and hence in the 
following 
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Y[(Mi-L + l) (2.3) 



will be used. /x(0,L) is the empirical probability density function (histogram) of local porosities. Its support is the 
unit interval. 

It is simple to determine fi((p , L) in the limits L — > and L — > oo of small and large measurement cells. For small 
cells one finds generally p5| , |f5[ 

/i(^,L = 0)=W-l) + (l-^) (2.4) 

where 

= V(PnS)/V(S) (2.5) 

is the bulk porosity. If the sample is macroscopically homogeneous then 

fjt(<f>, L -> oo) = 6(<t> - ~4>) (2.6) 

indicating that in both limits the geometrical information contained in fi(tfi, L) consists of the single number </>. The 
macroscopic limit, however, involves the question of macroscopic heterogeneity versus macroscopic homogeneity (for 
more information see flf|). In any case, if eqs. (2.4) and (2.6) hold it follows that there exists a special length scale 



L* defined as 

L* = min{L : /i(0, L) = /i(l, L) = 0} (2.7) 
at which the ^-distributions at 4> = and 0=1 both vanish for the first time. 

B. Local Percolation Probabilities 

The local percolation probabilities characterize the connectivity of measurement cells of a given local porosity. Let 

{1 : if K (x, L) percolates in "a" -direction 
(2.8) 
: otherwise 

be an indicator for percolation. What is meant by "a" -direction is summarized in Table I. 
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TABLE I: Legend for index a of local percolation probabilities X a (<j>,L). 



index a 


meaning 


X 


^-direction 


V 


y-direction 


z 


z-direction 


3 


(x A y A z)-direction 


c 


(xVyV z)-direction 





(-i(x V y V z))-direction 



A cell K(x, L) is called "percolating in the x-direction" if there exists a path inside the set P PI K(x, L) connecting 
those two faces of S that are vertical to the x-axis. Similarly for the other directions. Thus A3 = 1 indicates that the 
cell can be traversed along all 3 directions, while A c = 1 indicates that there exists at least one direction along which 
the block is percolating. Ao = 1 indicates a blocking cell. 

The local percolation probability in the "a" -direction is now defined through 

A a (^£)= Ex ^ (X ; £)W ' L) - ( 2 -9) 

2^x d ^0(x,L) 

The local percolation probability \ a (<fi, L) gives the fraction of measurement cells of sidelength L with local porosity 
4> that are percolating in the "a" -direction. 



C. Total Fraction of Percolating Cells 



The total fraction of all cells percolating along the "a" -direction is given by integration over all local porosities as 

Pa (L)= [ ^{<P,L)\ a {^L)d<P (2.10) 
Jo 

This quantitiy provides an important characteristic for network models. For a network model it gives the fraction of 
network elements (bond, sites etc.) which have to be permeable. 



III. SAMPLES AND ALGORITHMS 



The data sets of three different sandstones are used in the analysis below. Each data set consists of a threedimen- 
sional array of 0's and l's indicating pore space P or matrix M. The array dimensions are M\^Mi and M3. The 
pore space P of the three samples are displayed in Figures [l], || and ||. Note that the representations are not to scale 
because the resolution of each image is different. The data were obtained by computerized microtomography fl32| ]. 

Table II gives a synopsis of the characteristics of the three samples which have been analysed. 



TABLE II: Overview over properties of the data sets for three reservoir sandstones. 



Sample 


Description 


a 


Mi x M 2 x M 3 


L* 


4> 


k 


A 


Berea 


lOfim 


128 x 128 x 128 


260/mi 


0.1775 


HOOmD 


B 


Brent 


2.7/im 


180 x 217 x 217 


108/mi 


0.1602 


470mD 


C 


Sst20d 


30/im 


73 x 128 x 128 


540/im 


0.2470 


20000mD 
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Here a is the re solu tion, and Mi are the dimens ionl ess sidelengths of the sample in units of a. The bulk porosity was 
defined in eq. (2J3) and the length L* in eq. (2/7). The permeability k, given in millidarcy, is the experimentally 
determined permeability of the sample from which the data sets were obtained. 



The calculation of /i(0, L) is straightforward, and proceeds exactly according to eq. (2.2). Several possibilities exist 
for the choice of x in K(x, L). Originally |2^] it was proposed to choose for x a cubic lattice with lattice constant L 
such that 1J X K(x, L) = § and such that the resulting set of K(x, L) are nonoverlapping i.e. K(x, L) n K(x', L) = for 
x / x'. For the small data sets available this leads to poor statistics with strong fluctuations in all results. Therefore 
we use here a cubic lattice with smaller lattice constants giving rise to overlapping cells. The results below were 
obtained by using unit lattice constant. In other words we used for x all lattice sites except those whose distance 
from the sample boundary is less than L/2. It must be noted, however, that this method of positioning the cells 
gives progressively higher weight to the central region of the sample. This can, for large L, lead to small differences 
(roughly 0.005 in the present case) between the bulk porosity <j> as defined in eq. ( |2.5| ) and expected local porosity 
defined as f (j>fi(<j))d<j). 



The determination of A Q (x,L), i.e. of whether or not a cell is percolating in a given direction, was carried out 
according to the well known Hoshen-Kopelm an a lgorithm |33f] . The choice of x was the same as in the measurement 
of [i. \ a ((f>,L) was then calculated from eq. (2.9) 



IV. RESULTS 



The first sample is Berea sandstone whose pore space is displayed in Figure |l[ In this case the resolution is a = 10/xm, 
and the sidelengths of the sample are Mi = Mi = M3 = 128. Figure || shows the local porosity distributions /j,(<j>, L) 
with L = 40, 80, 120, 300^m exhibiting the typical crossover between the limits L = and L = 00. The curves are 
shown as dotted lines and marked by four different symbols corresponding to the four values of L as indicated in the 
legend. Next, in the same Figure^, the local percolation probabilities L) are displayed for the same four values 
of L that were used for fi. The curves for As(0, L) are distinguished from those for fj,(<j>, L) by a solid line style. The 
symbols used to indicate L are the same in both cases. The ordinate for the A3-graphs is the right axis, those for 
/i-graphs is the left axis. The local percolation probabilities A3 are increasing from zero to one. This expresses the 
fact that the full sample is connected. For L = 40/im= 4a the sidelength of the measurement cell corresponds to four 
voxels. To have a conducting path in all three directions one needs at least 10 voxels of pore space. This amounts 
to a porosity of roughly 0.16, and hence all curves A3(0,4a) must vanish below <j> « 0.16. Similarly to disconnect at 
least one of the three directions one needs at least 16 voxels filling a plane. Hence the curves A3 (0,4a) must equal 
unity above cf) « 1 — 0.25 = 0.75. This can be observed in Figure [| for Berea and in Figures || and [ll] for the other 
samples. For general L and dimension d the same consideration gives that A3 vanishes below (dL — d + 1)/L d and 
equals unity above 1 — 1/L. 



It is instructive to compare [i with A at a fixed L by superposing them in the same plot. Such a plot is shown 
in Figure || for L = L*. The characteristic length L* was defined in (2.7). For sample A (Berea) its value is found 
to be L* = 260/im. To facilitate comparison the local porosity distribution fj,(<fr,L*) has been rescaled such that its 
maximum equals unity. All six local percolation functions A Q are displayed in Figure ^. The sample appears to be 
isotropic because the three functions \ x ,\ y , X z all fall on top of each other. The curves A c and A3 are upper and lower 
bounds for the region inside which the connectivity increases from blocking to fully connected. Note also that this 
band is shifted to the left of the maximum of fi indicating that Berea sandstone is well connected. In fact the figure 
shows that an average cell (i.e. a cell with local porosity around 0.18) is percolating with probability larger than 0.75. 
Only cells with local porosity much below average are blocking. 



To investigate how heterogeneities in the connectivity are reflected in the local percolation probabilities, we have 
constructed an artificial modification of the Berea sample. To this end we have blocked roughly 1200 additional voxels 
out of the total of 128 3 w 2 • 10 6 voxels by blocking a plane of 1 x 100 x 100 voxels inside the sample. The orientation of 
the blocked plane was chosen perpendicular to the ^-direction, and the plane was centered in the y- and z-directions. 
It was placed in the middel, i.e. into the 64th layer along the x-direction. This amounts to a small decrease of porosity 
by roughly 0.00058. Note that the plane does not block the x-connectivity completely, but leaves an open shell at the 
sample boundary. The resulting modified porous microstructure was visually indistinguishable from the unmodified 
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one from the perspective of Figure [j]. Only when viewing the sample at right angles from the y- or z-direction it was 
possible to detect a small modification. In Figure || we display the same superposition of fi and A for the partially 
modified sample that was shown in Figure || for the unmodified sample. The local porosity distribution of the partially 
x-blocked sample is almost identical to that of the unmodified sample. This remains true for all values of L, and 
L* = 260/im is also unchanged. The functions Ao,A c ,Aj, and A z also remain almost unchanged. Aa; however differs 
from X y and X z as expected. As a consequence also A3 falls significantly below the result of the unmodified sample. 
The deviations of roughly 15% give an order of magnitude for the influence of connectivity fluctuations on A3. The 
difference between the results grows with increasing L. This allows another important conclusion. To characterize 
heterogeneities using local porosity analysis it is necessary to measure both [i and A as functions of L over a sufficiently 
wide range of L. Choosing only one fixed L may be misleading. 

The dependence of the difference between the unmodified Berea and the partially blocked sample on L has been 
f urthe r quantified in Figure 0. This figure shows the total fraction of percolating cells determined according to eq. 
( [2. 10 ). The results for the original unmodified sample are shown with solid lines, those for the sample with a partially 



blocking plane arc shown as dotted lines. This plot shows again that the unmodified sample is very isotropic because 
Px,Py,Pz overlap. In the partially blocked sample, however, deviations start to appear around L = 150/im in p x {L) 
and Pd,{L) becoming more pronounced at higher L. Note that the modified sample shows a decrease in p x {L) at 
L w 300/im while all curves for the unmodified sample are monotonously increasing. Of course p x (L) must start 
to increase again at L > 400/im because the sample is still connected on large scales. Therefore one expects that 
nonmonotonous behaviour of p(L) correlates with the length scale of heterogeneities in the connectivity. 

The next sample is Brent sandstone, shown in Figure ^ with a resolution of a = 2.7/im, and sample dimensions 
Mi = 180, M 2 = 217, M3 = 217. Although the data set for this sample is the largest one with respect to the number 
of voxels, its absolute size is the smallest of all samples. As a consequence the statistics of this sample is poor because 
it represents little more than a few pores. A larger sample seems necessary to obtain a representative sampling of the 
pore space. 

The local porosity distribution and local percolation probabilities for this sample are shown in Figure |8| using the 
same method of plotting as in Figure [| (see above). The superposition of fi((f>, L*) and \(<f>, L*) is displayed in Figure 
|j| Finally the total fraction of percolating cells is shown in Figure [H]. This curve seems to indicate that, while the 
sample is isotropic for small L, it shows increasing anisotropy at larger L. This effect may, however, also be due to 
the poor statistics as a result of the small absolute size of the system. 

Sample C is a clean weakly consolidated sandstone of unkown origin denoted below as Sst20d because its permeability 
is 20 Darcy. The sample has resolution a — 30/im and sample dimensions M\ = 73, M2 = 128, M3 = 128, and its pore 
space is displayed in Figure |^. The local porosity distribution for this sample is shown in Figure [ll] together with 
local percolation probabilities. Their superposition for L = L* is displayed in Figure |l^. Finally the total fraction of 
percolating cells is shown in Figure 

Figures |lj and [ll] indicate that the sample is anisotropic because A^ is significantly smaller than X y and X z . It has 
been checked that this is not a finite size effect due to the sidelength in the x-direction being shorter in this sample 
than, say, in sample A. This check was carried out by first dividing sample A in half along the x-direction, and then 
carrying out the analysis on the remaining truncated sample. This did not produce differences between X x , X y and X z 
for sample A. 

Hence it must be concluded that sample C is anisotropic in its connectivity, being less permeable in the x-direction 
than in the y- and z-directions. We point out, however, that to the unaided eye Figure || appears visually isotropic. 



V. COMPARISON OF RESULTS AND DISCUSSION 



Having presented the results for the various samples we now compare the samples against each other. Figure |l4| 
shows the local porosity distribution of all three samples at the same length L = 120/im. An important reason for 
the differences are different characteristic length scales for different samples. Sample C has clearly the largest length 
scale because its /i(0)-curve is closest to the L = limit of eq. (2.4). Next comes sample A, and sample B has the 
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smallest pores. The difference in characteristic length scales may be eliminated by comparing the samples at some 
intrinsic length scale such as the correlation length or the intrinsic length L* . 



Figure |l5] shows the comparison of all three samples at the intrinsic length L = L* . Now the local porosity 
distributions resemble each other much more closely. Nevertheless characteristic differences remain not only in their 
peak position but also in their shape. These may in part, but not entirely, be attributed to the different porosities. 
Sample A and sample B have nearly the same porosities, but the shape of their //(</>, L*) differs significantly. The width 
of the curves indicates the strength of porosity fluctuations, and hence is a quantitative measure of heterogeneities 
in the porosity. Using the width as a criterion we find that sample A is most homogeneous while sample B is 
most heterogeneous, and sample C is intermediate. This agrees with visual inspection of Figures [l] through and 
illustrates that n(4>) measures porosity heterogeneities. This fact was first demonstrated for twodimensional images in 
p3| and suggests to use local porosity distributions and percolation probabilities as a quantitative method to compare 
microstructures between models and experiment. 



In all three samples the behaviour of L) as function of L seems to approach the limits given in eqs. (2.4) and 
(2.6). This indicates that the samples approach macroscopic homogeneity for L —> oo |lq ]. Of course much larger 



samples (particularly for sample B) are needed to conclude this with certainty. 

The behaviour of X(4>, L) reflects the same trend towards macroscopic homogeneity because these functions approach 
a unit step at <f> with increasing L. The universality in the limit L — > is reflected in the fact that the A-curves for 
small L are very similar in Figures |, H and [tl]. This was already discussed above for L = 4a. Of course,for fixed L, 
the range of <j> over which A changes coincides with the range of 4> where n(4>, L) differs from zero. Figure fLQ compares 
the local percolation probabilities A3 and Ao for all samples at L = 120/im. Not surprisingly the three samples exhibit 
very different behaviour analogous to the difference in /i seen in Figure [l4|. 

Figure [l7] shows A3 and Ao of all samples at the intrinsic length scale L = L*. Characteristic differences in shape 
appear which emphasize the different connectivity of the three samples. These differences are not mere fluctuations 
because they are of the same order of magnitude as the differences introduced into A3 by the introduction of a blocking 
plane into sample A that was discussed above. The curves are statistically most reliable for A close to unity. This can 
be seen from Figures |[ ^| and ^|. These plots show that the maximum of /1 occurs at </>- values for which A is close to 
unity. They also show that the five functions X a with a = c, x, y, z, 3 fall to the left of the maximum of \x for sample A 
and C, while for sample B these functions change most rapidly in the vicinity of the maximum of fi. Thus for sample 
A and C the change from blocking to percolating occurs well below the average porosity in the low porosity tail of /1, 
indicating that both samples have a very high degree of connectivity. The modified sample A (with partially blocking 
plane) shows a broader overlap between Aq, and /i indicating lower connectivity (see Figure |^) , and sample B appears 
to have the lowest degree of connectivity. In summary Figure [l^ shows that although all the samples are very well 
connected the fluctuations in connectivity are different, and hence one must also expect permeability fluctuations. 

When comparing A3 for all samples it is seen that this function reaches a plateau at large <p. For sample A there 
is a region around <fi ~ 0.3 where A3 decreases (see Fig ||). This shows that local percolation probabilities are not 
always strictly monotonous (as might have been expected) but may exhibit minima and maxima indicating a variable 
fraction of blocking cells and hence connectivity heterogeneities at intermediate scales. 

Finally it is instructive to compare the total fraction of percolating cells p(L) for all samples. Figures [?|, [h] and 
O show that samples A and possibly B are essentially isotropic while sample C has clearly anisotropic connectivity. 



Figure [18] shows ps{L) as calculated from eq. (2.10) for all samples. This plot can be used for constructing network 
models in two ways. Firstly, if each K(x, L) is used to represent a site in the network model, p^,{L) provides an 
estimate for the fraction of percolating network units as a function of the networks lattice constant. Secondly, ps{L) 
exhibits intrinsic length scales. Extrapolating a tangent at the inflection point of each curve to p = 1 gives a length 
scale which could be interpreted as the minimum length scale for a representative elementary volume (REV) needed 
in homogenization and other averaging procedures ]3~i|Ji"5| ]. In this way one finds approximately 350/im for sample A, 
230/.tm for sample B and 700/mi for sample C as the smallest sidelength for an REV. 

The results presented here suggest that none of the three samples can be adequately modelled by a critical site 
percolation network. The reason is that in such a model the fraction of percolating sites should equal p% and hence 
would be strongly L-dependent, except for very small or very large L. Hence, plots such as Figures ^ or provide 
information on how to choose the network elements (site, bonds) of a network model and how to relate the length 
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scale of the real rock to the networks lattice constant. 

Figure |l^ shows ps(L/L*) for all samples as solid lines. The dotted line corresponds to the modified sample A with 
a partially blocking plane. The samples are again significantly different. The vertical shift is in part due to differences 
in porosity because generally ps{L = 0) = <f>. The sample with the blocking plane indicates that it is possible to 
have intermediate plateaus and nonmonotonicity, at least if large scale heterogeneities are present. The unmodified 
homogeneous samples exhibit different widths over which the curves increase from ip to unity. These widths may be 
used as a quantitative measure of fluctuations in connectivity, both in theoretical models and in experiment. 
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FIGURES 



FIG. 1. Threedimensional pore space of Berea sandstone (sample A). The resolution is a — 10/im, the sample dimensions 
are Mi = 128, M2 = 128, M3 = 128. The bulk porosity is <j> — 0.1775. The pore space is indicated in blue, the matrix space is 
transparent. 

FIG. 2. Threedimensional pore space reconstruction of Brent sandstone (sample B). The resolution is a = 2.7/xm, the sample 
dimensions are Mi = 180, M2 = 217, M3 = 217. The bulk porosity is 4> — 0.1602. The pore space is indicated in blue, the 
matrix space is transparent. 

FIG. 3. Threedimensional pore space reconstruction of a weakly consolidated sandstone (sample C). The resolution is 
a — 30/im, the sample dimensions are Mi = 73, M2 = 128, M3 = 128. The bulk porosity is cj> — 0.2450. The pore 
space is indicated in blue, the matrix space is transparent. 
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FIG. 4. Local porosity distributions (j,(cj>,L) (dotted lines) and local percolation probabilities \^{4>,L) (solid lines) for Berea 
sandstone (sample A, cf. Figure [l]). Four different values for L are indicated by different symbols defined in the legend. The 
ordinate for the graphs of jj,((j>,L) is on the left, the ordinate for \z{(f),L) is on the right as indicated by the axis labels. 
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FIG. 5. Local percolation probabilities \ a {4>,L*) with a — 0, 3, c, x, y, z and L* — 260^tm for Berea sandstone (sample A) 
shown in Figure |fj. The dotted curve with circular symbols is the local porosity distribution at L = L* — 260/xm rescaled to 
have maximum 1. 
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FIG. 6. Same as Figure but for sample A modified with partially blocking plane in i-direction. See text for details. 
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FIG. 7. Total fraction of percolating cells p a (L) with a — 0,3, c, x, y, z for Berea sandstone (sample A, cf. Figure [I]) shown 
as solid lines, and for the modified sample A (with partially blocking yz-plane) shown as dotted lines. 




FIG. 8. Local porosity distributions fj,((f>,L) (dotted lines) and local percolation probabilities \2,{4>,L) (solid lines) for Brent 
sandstone (sample B) shown in Figure ^| Four different values for L are indicated by different symbols defined in the legend. 
The ordinate for the graphs of n{4>,L) is on the left, the ordinate for \^{4>,L) is on the right as indicated by the axis labels. 
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FIG. 9. Local percolation probabilities \ a (<f>,L*) with a — 0, 3, c, x, y, z and L* = 108/im for Brent sandstone (sample B) 
shown in Figure g. The dotted curve with circular symbols is the local porosity distribution at L — L* = 108/im rescaled to 
have maximum 1. 
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FIG. 11. Local porosity distributions _u(<f>,L) (dotted lines) and local percolation probabilities \z(<f>,L) (solid lines) for the 
sandstone (sample C) shown in Figure pi Four different values for L are indicated by different symbols defined in the legend. 
The ordinate for the graphs of fj,(<f>,L) is on the left, the ordinate for A3 (<j>,L) is on the right as indicated by the axis labels. 
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FIG. 15. Local porosity distributions n(4>,L — V) for all three samples 
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FIG. 16. Local percolation probabilities \z(<f),L — 120/im) and Xo(4>,L = 120/xm) for all three samples 
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FIG. 17. Local percolation probabilities \-s(<f),L = L*) and \o(<j>,L — L*) for all three samples 
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